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Abstract. We investigate the possibility of realizing unconventional superconductiv- 
ity in doped band insulators on the square lattice and the honeycomb lattice, where 
the latter is found to be a good candidate due to the disconnectivity of the Fermi 
surface. We propose applying the theory to the superconductivity in doped layered 
nitride /3-MNCl (M =Hf, Zr). Finally, we compare two groups of superconductors 
with disconnected Fermi surface, /3-MNCl and the iron pnictides, which have high 
despite some faults against superconductivity are present. 
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1. Introduction 

Layered nitride /3-MNCl[T] (M =Hf, Zr) doped with carriers is one of the most 
interesting group of superconductors. The mother compound /3-MNCl is composed of 
alternate stacking of honeycomb MN bilayer and CI bilayer. [2] This is a band insulator 
and becomes a superconductor upon doping electrons by Na or Li intercalation. They 
exhibit relatively high Tc up to ~ 25K for M =Hf and ~ 15K for M =Zr. The bilayer 
honeycomb lattice structure consisting of M and N (Fig{T]) is considered to be playing 
the main role in the occurrence of superconductivity, and the two dimensional nature of 
the superconductivity has been revealed by nuclear magnetic resonance |3] and muon spin 
relaxation studies. [H [5] Despite the relatively high Tc, experimental as well as theoretical 
studies indicate extremely low density of states (DOS) at the Fermi level. P, [71 E] In 
fact, they have the highest among materials with the speciflc heat 7 as small as ~ 
ImJ/molK^. The electron phonon coupling is also estimated to be weak,[6l El E], [ID] 
and the isotope effect is found to be small. [HI |12] In the superconducting state, the 
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density of states recover rapidly upon increasing the magnetic field, [7] suggesting some 
kind of anisotropic pairing. As for the doping dependence, the DOS at the Fermi level 
stays nearly constant, but, for Lia;ZrNCl, Tc shows an increase upon lowering the carrier 
concentration until a sudden disappearance of the and a superconductor-insulator 
transition is observed. |13] On the other hand, in Lia;HfNCl, Tc stays nearly constant 
for the whole doping range x < 0.5. [Hj Furthermore, for Li^HfNCl, an intercalation of 
organic molecules tetrahydrofuran (THF) between the layers is found to enhance T^. [H] 
While these experiments suggest that some kind of unconventional pairing mechanism 
may be at work, tunneling spectroscopy experiments on the other hand find an s-wave 
like, fully open gap. [15] Specific heat measurements also suggest an s-wave like gap, but 
again the doping dependence of the magnitude of the gap is unusual. In the underdoped 
regime, the gap is large, while the gap becomes small as the doping level is increased, 
varying from a "strong coupling" to an "extremely weak coupling" superconductor. p!7] 




Figure 1. The bilayer honeycomb lattice is shown. The black and white circles 
represent M(=Hf, Zr) and N, respectively. 

In the present study, we first generally study the possibility of unconventional 
superconductivity by doping band insulators on a square or honeycomb lattice. [16] We 
find that the honeycomb lattice is a good candidate for realizing such a possibility, where 
the disconnect ivity of the Fermi surface plays a role. Secondly, we consider applying 
the theory to /3-MNCl. The two bands sitting closest to the Fermi level in the first 
principles band calculation [S] can roughly be reproduced by a single honeycomb lattice 
model, where the above general theory can be applied. Finally, we compare two groups 
of superconductors with disconnected Fermi surface, /3-MNCl and the iron pnictides, 
which have high despite some faults against superconductivity are present. 
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Figure 2. An explanation of tiie d-wave superconductivity due to spin fluctuations 
on the square lattice. The arrow represents the wave vector Q at which the spin 
fluctuations develop. The solid (dashed) lines are the portions of the Fermi surface 
where the gap has a positive (negative) sign. 



2. Square lattice 

Hubbard Model is a model that considers the on-site repulsive interaction f/ on a 
tight binding model. Let us start with the Hubbard model on a square lattice, where 
we consider only the nearest neighbor t as the hopping integral. When the band 
filling n (=number of electrons/number of sites) is near half-filling (n ~ 1), strong 
antiferromagnetic spin fiuctuations arise, and the possibility of d-wave superconductivity 
mediated by these spin fiuctuations has been discussed for the past several decades. The 
(i-wave superconductivity can be understood as follows. Generally, superconductivity 
occurs due to pair scattering mediated by the pairing interaction V{q). The gap equation 
can be written as 

Aw=-i: '"'''fj|;;y^^V (/.--fc)A(fe'). (1) 

where A is the gap, E is the quasiparticle dispersion, and T is the temperature. When 
the spin singlet pairing interaction is mediated by spin fiuctuations, V{k) is positive 
and takes large values at large Q, where the spin fiuctuation develops. So in the 
case of the square lattice, where the spin fiuctuations develop near Q = {n,7r), we 
have to change the sign of the gap between the wavevectors ~ (vr, 0) and ~ (0, tt) in 
order to have a finite A as a solution for the gap equation, which results in a (i-wave 
gap as shown in Figj2l By applying fiuctuation exchange (FLEX) approximation to 
this system, which is a kind of self-consistent perturbation theory that collects random 
phase approximation type diagrams, we can obtain the Green's function and the spin 
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susceptibility. [IB] These can be plugged into the Ehashberg equation, whose solution 
gives a d-wave superconductivity with a of the order of O.Olt, where t is the nearest 
neighbor hopping integral (if t ~ leV, Tc is of the order of lOOK). 

Now, the square lattice (with only the nearest neighbor hopping) is a bipartite 
lattice which can be separated into A and B sublattices. Let us see what happens 
to the superconductivity if we introduce a level offset A between A and B lattices. 
The introduction of A opens up a gap at the center of the band, so this amounts 
to investigating the possibility of unconventional superconductivity by doping carriers 
in band insulators (Fig|3](a)). As shown in FigJU we find that the introduction of A 
rapidly suppresses superconductivity. Thus, in the case of the square lattice, chances 
for realizing unconventional superconductivity in the above sense are small. 




Figure 3. (a) The density of states for the square lattice with A = (black) and 
A = t (red).(b) The density of states (upper) and the band dispersion (lower) for the 
honeycomb lattice with A = (black) and A — 1.5t. The hexagonal Brillouin zone of 
the honeycomb lattice is shown in the inset of the lower panel. 



3. Honeycomb lattice 

Now, let us compare the above result for the square lattice with those for another two 
dimensional bipartite lattice, namely, the honeycomb lattice. The honeycomb lattice 
is a system with two sites in a unit cell, in which the two bands make point contact 
at K and K' points of the Brillouin zone, resulting in a zero gap density of states 
(Figini^b)). We show in Figin](a) the maximum value of the spin susceptibihty as a 
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Figure 4. Tc as functions of the level offset A obtained by FLEX+Eliashberg equation 
for the square lattice with U — 6t, n ^ 0.7 (black) or for the honeycomb lattice with 
U = 6t and n = 1.08 (red). 




Figure 5. (a) The maximum value of the spin susceptibility as functions of 
temperature obtained by FLEX for the square lattice with n = 0.7 or n = 0.65 and 
for the honeycomb lattice for n ~ 1.08. U = 6t in all cases. HC and SQ stand for the 
honeycomb and the square lattices, respectively, (b) The eigenvalue of the linearized 
Eliashberg equation with the same parameter values as in (a). The temperature at 
which A = 1 is the T^. 



function of temperature for the band filling oi n = 1.08, which corresponds to a small 
amount of electron doping. Surprisingly, we find that the spin susceptibility is nearly 
independent of T, which is in sharp contrast with the case of the square lattice. For 
example, for the square lattice with n = 0.7, which is already substantially away from 
half-filling, we still have a strong enhancement of the spin susceptibility upon lowering 
the temperature. With further hole doping to n = 0.65, the spin susceptibility is 
suppressed, but even in that case, there is a moderate increase of the spin susceptibihty 
upon lowering the temperature. In Figl5]^b), we show the eigenvalue A of the linearized 
Eliashberg equation as functions of temperature. Tc is the temperature where A reaches 
unity. The density of states at Ep is nearly the same for the square lattice with n = 0.65 
and the honeycomb lattice with n = 1.08, and also the spin susceptibility has similar 
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Figure 6. The contour plots of the FLEX resuh at the lowest Matsubara frequency 
for the honeycomb lattice in the hexagonal Brillouin zone with U — 6t, n — 1.08, 
and T — O.Oli (a) The Green's function of the upper band squared, (b) the spin 
susceptibility, (c) the superconducting gap function. 




Figure 7. The Fermi surface (the two circles) and the sign of the gap function are 
schematically shown in the extended zone scheme. The dashed arrows represent the 
wave vectors at which the spin fluctuations develop. 

values at low temperature, but still, the honeycomb lattice has larger A and higher Tc. 
Thus, the Hubbard model on the honeycomb lattice has relatively high despite the 
low density of states and the weak and temperature independent spin fluctuations. 

Figj6] shows the contour plot of the Green's function squared, whose ridge 
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corresponds to the Fermi surface. We see here two disconnected pieces of the Fermi 
surface. The spin susceptibihty is maximized at wave vectors that bridge the opposite 
sides of each pieces of the Fermi surface. As can be seen more clearly in FigJTl the 
gap has a d-wave form, where the gap changes sign across the wave vector at which 
the spin susceptibility is maximized. Here, note that one of the nodes of the gap do 
not intersect the Fermi surface because of its disconnectivity, which may be one reason 
why superconductivity is favored despite the low density of states and the weak spin 
fluctuations. By symmetry, there are two degenerate d-wave gaps (say, d^y and dj.2_y2, 
or any two linearly independent combinations), and the most probable form of the gap 
below Tc is the form d + id, where the two d-wave gaps mix with a phase shift of 7r/2. 
Since the two ci-wave gaps have nodal lines at different positions, this kind of mixture 
leads to a gap that has a finite absolute value on the entire Fermi surface. An important 
point is that if such a state is realized, the time reversal symmetry should be broken. 

Now we introduce the level offset between A and B sites as we did for the square 
lattice. In this case also, a band gap opens in the center (FigOJ^b)), so we once again 
investigate the possibility of superconductivity by doping band insulators. In this case, 
we find that superconductivity is relatively robust against the introduction of A. This 
may be because the density of states is already low in the original honeycomb lattice, 
so that the introduction of A does not affect superconductivity so much. 

4. Application to /?-MNCl 

We now consider applying the above theory to /3-MNCl. Although /3-MNCl has a bilayer 
honeycomb lattice structure, we find that the two bands closest to Ep obtained in the 
first principles calculation [8} [lOl [19], [20] can roughly be reproduced by a single layer 
honeycomb lattice model consisting of alternating "M" and "N" orbitals with a level 
offset as shown in FiglHl Here we take t = 1.2eV, A/t = 2.7, and t' /t = 0.35. If we 
consider on-site repulsive interaction [/ = 6t on both M and N orbitals, [21] the model 
is similar to the one studied in section 2, except that distant hopping integrals t' have 
to be considered so as to reproduce the first principles band structure. Consequently, 
within the FLEX+Eliashberg equation approach, we obtain relatively high of around 
30K as shown in Fig|9l 

Now we compare the present scenario with the experiments for /3-MNCl. As 
mentioned in the Introduction, relatively high is obtained despite the extremely 
low density of states and the weak electron phonon coupling, [6l [Tj, [91 [HI [ID] which 
can be explained within the present theory. The isotope effect is small, [TT| [T2] which 
again seems consistent since the the present pairing mechanism is purely electronic. 
The cointercalation of THF molecules enhances TcjH] and this also seems to be 
understandable within this kind of spin fluctuation mediated pairing, where the quasi 
two dimensionality is favored as discussed generally in refs. |22l 123] . As for the pairing 
symmetry, a fully open gap is observed in the experiments. |15j and this is consistent with 
the present scenario provided that the d + id state is realized. It is hence interesting to 




Figure 8. Upper panel: the tightbinding model considered for /3-MNCl. Lower panel: 
the first principles band calculation taken from ref . [8] and the band dispersion of the 
tightbinding model. 

investigate experimentally the possibility of time reversal symmetry breaking in the 
superconducting state of this material. The rapid recovery of the specific heat by 
applying the magnetic field, [7] and also the unusual doping dependence of both the 
Tc and the magnitude of the gap[T3l [HI [T7j remain as an interesting future problem. 
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Figure 9. Tc plotted as a function of the band filling for the model shown in Fig[Hl 

5. Superconductivity in systems with disconnected Fermi surface 

Finally, let us go back to the d-wave superconductivity on the square lattice. As 
mentioned, Tc obtained by FLEX is the order of O.Olt, where t is the hopping integral. 
This can correspond to a high in actual materials because t can be of the order of 
electron volt, but still this is low compared to the original energy scale t. There are 
several reasons for the "low" T^, and one of them is that we have to have nodes of the 
gap intersecting the Fermi surface because sign change of the gap is required in the case 
of spin fluctuation mediated pairing. In this context, we proposed some time ago that 
if there are disconnected pieces of Fermi surface nested to some extent, we can change 
the sign of the gap between the disconnected peaces without the nodes intersecting the 
Fermi surface, and this can result in a high Tc superconductivity. |24j An example of 
such a Fermi surface is shown in FigJTOl Possibilities of the disconnected Fermi surface 
playing important role in the occurrence of superconductivity or the determination of 
the pairing symmetry have been discussed for a cobaltate Naa;.Co02[25] and an organic 
superconductor (TMTSF)2X.[26] 

Quite recently, superconductivity has been found in iron-based pnictides by Hosono 
and coworkers. [27] The highest Tc of this series of materials have reached 55K.[28] Band 
calculations show that there are several disconnected pieces of the Fermi surface in 
this material. |29j where spin fluctuations can arise due to the nesting between them. [30] 
According to our Eliashberg theory calculation that takes into account such kind of spin 
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Figure 10. An example of fully open sign reversing gap on a disconnected Fermi 
surface is shown. The arrow shows the nesting vector at which the spin fluctuations 
develop. 

fluctuations, [Slj the gap changes sign across the nesting vector of the Fermi surface, and 
the magnitude of the gap is especially large on the portion of the Fermi surface where 
the dx2-y2 orbital character is strong as shown in FigJTTl The nearest neighbor hopping 
of this orbital is about O.lSeV, which is quite small, and the experimentally observed 
maximum = 55K corresponds to about 0.04t which is higher than that can be reached 
in the single band square lattice. 

So the two groups of materials seem to exhibit high despite some kind of faults 
against superconductivity are present: for /3-MNCl, relatively high Tc is obtained despite 
the extremely low DOS and the weak spin fluctuations (and electron-phonon coupling), 
while for the pnictides, high Tc is obtained despite the low energy scale of the main 
band, and also competing spin and charge fluctuations due to the multiplicity of the 
orbitals. The disconnectivity of the Fermi surface may be one good reason why these 
faults are overcome. In the future, there may a possibility that we can get higher by 
realizing disconnected Fermi surface on more ideal situations. 

The author acknowledges Y. Taguchi, Y. Iwasa, Y. Kasahara, T. Ekino, and H. 
Aoki for fruitful discussions. Numerical calculation has been done at the computer 
center, ISSP, University of Tokyo. This study has been supported by Grants-in Aid 
from MEXT of Japan and JSPS. 
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